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We point out that double distributions need not vanish at their boundary. Boundary terms do 
not change the ambiguity inherent in defining double distributions; instead, boundary conditions 
must be satisfied in order to switch between different decompositions. We analyze both the spin 
zero and spin one-half cases. 
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QCD factorization provides the way to access infor- 
mation experimentally about the non-perturbative quark 
and gluon substructure of hadrons. In recent years, much 
attention has been generated by hard exclusive reactions 
such as deeply virtual Compton scattering, and the hard 
electroproduction of mesons The non-perturbative 
structure functions entering these reactions are general- 
ized parton distributions (GPDs), see the reviews @- 

The polynomiality property required of the Mellin mo- 
ments of GPDs is elegantly explained by the formalism of 
double distributions (DDs) [3J, |j, |5j . Phenomenological 
modeling of GPDs is almost exclusively done utilizing 
parametrized DDs. Experimentally one cannot access 
the DDs directly, only indirectly through the H(x,^,t) 
and E(x, £, t) GPDs. In Ref. [||, the ambiguity inherent 
in defining DDs for the pion was likened to the gauge 
ambiguity of the vector potential of a two-dimensional 
magnetic field. Here we extend this analysis to the pro- 
ton case. But first, we review the pion case and tie up 
the loose ends relating to the non- vanishing value of the 
DDs at their boundary. Additionally we compare various 
forms of DDs. 

For the pion, we have the following decomposition of 
the non-diagonal twist-two matrix elements in terms of 
various form factors A n k(t) and B n k(t) 

(P' ^(O^iD^ ■ ■ ■ iD^ } ijj(0)\P) 

xpw -^(-i) -(-2) ■ (1) 

where the gauge covariant derivative D = (D — D)/2, 
{. . .} denotes the symmetrization and trace subtraction 
performed on Lorentz indices, P = (P' + P)/2, A = 
P' — P, and t = A 2 . Often we treat the t-dependence as 
implicit below. Above, T-invariance restricts A n k(t) = 
for k odd, and B n k(t) = for k even. There is mani- 
fest arbitrarity in the twist-two form factors appearing 
in Eq. 0). The particular decomposition above can be 
used to define two DDs for the pion. These DDs are 



generating functions for the twist-two form factors 

Above we have abbreviated the integration as a = 

d(3 J^jjL da. As a consequence of T-invariance, the 
function F([3, a) is even in a, while G(/3, a) is odd. 

Summing up the moments in Eq. ifTJl. these DD func- 
tions then appear in matrix elements of the light-like 
separated quark bilinear operator A4(P ■ z, A • z) = 
(P'\Jj(-z/2)^(z/2)\P), 

M(P-z,A-z)= [ e -tfF-*+*«A-»/2 

x [2P ■ zF{P,a) - A- zG((3,a)] , (3) 
where z 2 — 0. Now we define the pion GPD 

H(x,0 = lJ^e ix7+z -M(P-z,A.z), (4) 

with the usual definition £ = — A+/2P+. 

Any other choice of generating functions for the mo- 
ments must lead to the same GPD. To expose this am- 
biguity, we follow 0. Integration of Eq. J2J by parts 
produces surface terms that in general do not vanish. 
We are careful about this point because previous mod- 
els 0] indeed have non- vanishing contributions on the 
boundary |/3| = I — Mathematically the DDs must 
vanish only at the corners of support: <5(a)<J(|/3| — I), 
<5(/3)<$(|a| — I), else the form factors in Eq. J2J) do not 
fall off in the space of moments. Physically the DDs' 
vanishing at the first set of corners is tied via Eq. Q to 
the vanishing of the GPDs at x = ±1, which is known 
from perturbative QCD . Vanishing of the DDs at the 
second set of corners implies the continuity of GPDs at 
the crossover (x — £), which in turn ensures factoriza- 
tion. Boundary contributions are thus not ruled out and 
the DDs appear to be loose at the ends. Phenomeno- 
logical parameterizations of DDs must generally include 
such terms. The complete result of integrating by parts 
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can be expressed in the form 

M(F-z,A-z) = -2z [ N(p,a)e- l/3T! - z+laA - z ^ 

J0,a 

+8i f da cos(aA ■ z/2) 
Jo 

x jcos [(1 - ajP-z] S + (a) 
-isin[(l-a)P-z] S~{a)\, (5) 

where N((3,a) = -§pF{f3,a) + ^G(/3,a) and S ± (a) = 
F+(l — a, a) + G ± (l — a, a). In deriving Eq. J3J), we made 
use of the a-symmetry of DDs and additionally have used 
i-distributions (F ± and G ± ), which are even and odd 
functions of j3 1 respectively. 

Consider an arbitrary potential function x(Pi a ) which 
is odd with respect to a. We can decompose x m terms of 
its even and odd parts with respect to /3, namely x^ {P, 
where x^(~ Pi a ) — ^X^iPi®)- The potential function 
can be used to generate the DD transformation 



G ± (P,a) 



G±(A a )-|f(Aa)j- W d 



Notice this transformation preserves both the a- and 
/3-symmetry of the i-distributions. Next we observe 
Eq. (jSJ is invariant under this transformation provided 
the potential also satisfies two boundary conditions 



d_ 

da 



0=1- 



(7) 



0=i- 



Above both a and f3 are positive and the corresponding 
boundary conditions for the full range of DD variables 
can be found trivially due to the symmetry properties. 
Without the boundary conditions Eq. J7|) on the potential 
function, the transformation generated by Eq. © is that 
of 0. 

To convert the DDs in Eq. © to the Polyakov- Weiss 
"gauge" Q , in which there is only an F-type double dis- 
tribution and D-term, we use the potential specified by 



xf(/J.«) = -o 



i+M 



dp' G^{P',a) 



1 ° dp' G^(P', a) - I (1 T 1) sign(P)D(a) 



where D(a) is the _D-term given by 

rl-\u\ 



(8) 



D(a) 



dpG+(P,a). 



(9) 



One can verify that the potential specified by Eq. © 
satisfies the boundary conditions Eq. Q provided one 



also assumes the GPD is continuous, i.e., D(±l;t) = 
0. Under the transformation generated by Eq. (JSJ), the 
bilocal matrix element reads 



M(P-z,A-z) = e -i0P-z+iaA-z/2 
J0,a 

x [2P ■ z F (p,a) - A- z S(p)D(a)] , (10) 
where the resulting DD F a (P, a; t) is given by 

F ((3, a) = F(p, a) + ^ [ X + {P, a) + X ' {P, a)} . (11) 

A perhaps more interesting choice of gauge is what we 
call the Drell-Yan gauge. It is specified by the potential 
function 



xt(P,a) = - 1 - 

1-1/31 



1+1/31 



da'F ± (P,a') 



da' F (P,a') - sign(a) 7J (/3) 



, (12) 



where the Z?-term D(P) is given by D(P) — D + (P) + 
D-(P), with 



D ± (P) 



1-1/31 



1 + 1/31 



daF ± (P,a). 



(13) 



Again one can verify that the potential specified by 
Eq. l(T2l satisfies the boundary conditions Eq. Q) . Under 
this transformation, the bilocal matrix element reads 



M(P-Z,A- Z ) = e -i0P-z+i<*A-z/2 
J0,a 

x [2P-z 6(a) D(P) - A-z G x (/?,«)] , (14) 
where the resulting DD Gi(/3, a) is given by 

Gr(/3, a) = G(J3, a) - ^- [*+(/?, a) + X i(P, <*)] • (15) 

The Drell-Yan gauge is particularly interesting from 
the perspective of GPDs. Inserting Eq. (|14(1 into the def- 
inition of the pion GPD Eq. (@J , we have (now reinstating 
the ^dependence) 



H(x,i,t) =D(x;t) +£ / 8(x- P-£a)Gi{P,a;t). 

J0,a 

(16) 

Notice the contribution to the GPD from the Z?-term is 
independent of £ and the reduction relations are simply 

h(x)=H(x,0,0) = D(x,0), (17) 

for the quark distribution and 

F(t)= dxH(x,£,t)= dxD(x;t), (18) 
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for the pion form factor. Thus the ZJ-term has a sim- 
ple physical interpretation as the x-integrand of the form 
factor calculated in the Drell-Yan frame (£ = 0). Chang- 
ing to a frame where £ ^ 0, we have Eq. Note, the 
second term in the GPD proportional to £ contributes 
nothing to the form factor because of Lorentz invariance. 
This is maintained by the a-symmetry of the DD. The 
GPD, by contrast, is manifestly affected by the change in 
frame. At fixed t, the change to £ ^ requires a dynam- 
ical light-front rotation. Despite claims ||, the D-term 
is all that one can learn about DDs from the Drell-Yan 
expression for the form factor |ic) . 

To describe the pion GPD, we have resorted to using 
two DDs F and G because these are encountered in ac- 
tual calculations Q. The GPD, however, can be viewed 
as the projection of a single DD function see also 
. To see this representation one must use a non-trivial 
transformation generated by X/(/3, a), where 



X/(0,a)=XoGM)+a / d/3'W(/3, /3')/ T (/3', «), 

J — l+\a\ 

(19) 

with W{P,(3') = 0{p)6{p - (3) - 8(-P)0(P - &). The 
new DD /(/?, a), which is even in a, is implicitly defined 
by the result of the transformation 



M(P-Z,A- Z ) = e -WP-z+iaA-z/2 
J/3,a 

x [2P-Z0- A-za] f(0,a). (20) 
The boundary conditions are met provided 

f*(l-a,a) = F?(l-a,a) ) (21) 

where F a is given in Eq. Ijlljl. 

Now we address the ambiguities of proton DDs. The 
non-diagonal proton matrix elements of twist-two opera- 
tors can be decomposed into form factors A n k (t) , B n k (t) 
and C nk (t) 

(P', A'^Ofr^iB" 1 • • ■ tB"»ty(0)|P, A) 



= «v(P')E 



fc=0 

ia^ v A 



k\{n-k)\ 



2M 



A^ 

— B n k{t) -—rG n k{t) 



AM 



MP) 



A 



(22) 



T-invariance forces A n k{t) = B n k(t) = for k odd, and 
C n k(t) = for k even. There are three Dirac structures 
in the above decomposition since in general the twist-two 
currents are not conserved. 

The above decomposition can be used to define three 
double distributions as generating functions for the twist- 



two form factors 




0,a 



~ k a k | K(/3,a) 



(23) 



T-invariance implies the functions P(/3, a) and K((3,a) 
are even in a, while G(f3, a) is odd. Eq. (1221) represents a 
"physical" gauge for proton DDs since these are encoun- 
tered in actual calculations [l2|. 

Summing up the moments in Eq. Q22[) . the DD func- 
tions then appear in matrix elements of the light-like 
separated quark bilinear operator A4 X ' X (P • z, A ■ z) = 
(P',X'\^(-z/2)^(z/2)\P,X), 



M X > X (P ■ z,A- z) 



-if3P-z+iaA-z/2— 



Uy(P') 



J,T?!R \ , «7 f "%A J/ A-z 
2M Ki ^ a) -Uf 



(24) 



Now we define the light-cone correlation function 



M x '' x (x,0 



dz -n 



2tt 



P+z- 



M x '> x (P-z,A-z). (25) 



This correlation function can be written in terms of the 
two independent GPDs H(x,£) and E(x,t;) 



M x ^ x (x,0 



1 



2P 



■ux'(P') 



u x (P). (26) 



Inserting the DD decomposition Eq. Q24JI into the corre- 
lator in Eq. 1)25(1 . we can express the GPDs as projections 
of the DDs 



H(x,0 
E(x,0 



(27) 

from which we can view the ^-dependence of GPDs as 
arising from different slices of Lorentz invariant DDs. 
Due to the symmetry of the DDs with respect to a, the 
GPDs H(x, £, t) and E(x, £, t) are both even functions of 
the skewness parameter £. 

Utilizing the Gordon identities, we can rewrite the bilo- 
cal matrix element in Eq. (|24|l as 



M x ' A (P-z,A-z) 



uy(P') 

2M ~ 2M J l3,a 



-if3P-z+iaA-z/2 



2P ■ zG E (l3,a) - A ■ zG(/3,a) 



i^ Zll A v P al ^G M {p,a) 



ux{P), (28) 
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where we have defined new double distributions 

G B (j3,a)=F((3,a) + ^—K(l3,a) (29) 
G M (l3,a)=F(P,a)+K(/3,a) (30) 
G(/3,a) = i(l-t/4M 2 )G(/3,a) (31) 

in analogy with the Sachs electric and magnetic form 
factors. After integrating by parts, we have 



M y ' X (P-z,A-z) 



iuy(P') ( f 



-i0P-z+iaA-z/2 



2N(f3, a) + e^z^A,P al0l5 G M (P, a) 



da cos 



(aA • z/2)|cos [(1 — a)P ■ z\ S + (a) 
i sin [(1 - a)P ■ z] S~ (a)} J u\(P), (32) 



where N((3,a) = &G E (P,a) + ^G(/3,a), and the 

boundary terms 5 ,± (a) = — a,a) + Cr (1 — a, a). 

Now we define a potential function x([3, a) as above. 
The expression Eq. (|32() is invariant under the transfor- 
mation 



G%{l3,a) 
,G±(/3,a) 



'G% 09, a) + 09, a) 

v G ± (/3,a)-^x T (/?,a) 



(33) 



provided the boundary conditions on x Eq. Q are sat- 
isfied. Translating the transformation Eq. Ij33(l into the 
original DDs, we have invariance under 



/ F(J3, a) 
\G(j3,a) 



\G(p,a)-&,x(P,<*) 



(34) 



As with the pion case, one can convert the proton DDs 
to Polyakov- Weiss gauge or Drell-Yan gauge using the 
potentials xtiPi a ) an d Xx(Pi a )i respectively. Addi- 
tionally two independent DDs for the proton can be un- 
masked using the transformation generated by X/O^j °)- 
Each of these transformations for the proton case is un- 
derstood under the replacement {F(0, a) — > Ge{$, £*)} 
and {G(/3, a) — > G(/3, a)}. In Polyakov- Weiss gauge we 
recover the result of [4| for the proton. In Drell-Yan 
gauge, we generate an additive ^-independent contribu- 
tion to the combination of GPDs H(x,£) + j^E{x,^)\ 
this is the x-integrand of the Sachs electric form factor in 
the Drell-Yan frame. The £ dependence is contained in a 
term analogous to the one in Eq. (|16|l and characterizes 
the dynamical light-front rotation to £ ^ 0. In the min- 
imal gauge generated by X/(A <*)■, we see there are only 
two underlying independent DDs. 



There are two additional Drell-Yan type gauges for the 
proton. Using Xi{j3, a) as appears in Eq. (JT^Jl, we gen- 
erate an additive ^-independent contribution to iJ(a;,£). 
This is the x-integrand of the Dirac form factor calculated 
in the Drell-Yan frame. Alternately we can use Eq. H12|) 
under the replacement {F(/3, a) — * —K(/3,a)} to gener- 
ate the £- independent, x-integrand of the Pauli form fac- 
tor in the Drell-Yan frame, which contributes to E(x, £). 
With Eq. Ij34(l . we cannot simultaneously find Drell-Yan 
contributions to GPDs from both Dirac and Pauli form 
factors. This limitation arises from Eq. (|33|) . as we can- 
not generate an additive contribution to H(x, £)+E(x, £) 
from the x-integrand of the Sachs magnetic form factor. 
Thus Eq. I|32l) may contain additional freedom. 

Above we have seen that DDs do not necessarily vanish 
at their boundary and have investigated the consequences 
of boundary terms on the ambiguity inherent to DDs. 
The potential function that generates a transformation 
of DDs must satisfy boundary conditions. We have car- 
ried out this analysis for both the spin zero and spin 
one-half cases. Inclusion of boundary terms into model 
DDs is necessary for general phenomenological parame- 
terizations of GPDs. 

We thank W. Detmold and G. A. Miller for useful com- 
ments. This work was funded by the U. S. Department 
of Energy, grant: DE-FG03 - 97ER41014. 
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